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We theoretically study the interaction of a heavy hole with nuclear spins in a quasi-two- 
dimensional III-V semiconductor quantum dot and the resulting dephasing of heavy-hole spin states. 
It has frequently been stated in the literature that heavy holes have a negligible interaction with 
nuclear spins. We show that this is not the case. In contrast, the interaction can be rather strong 
and will be the dominant source of decoherence in some cases. We also show that for unstrained 
quantum dots the form of the interaction is Ising-like, resulting in unique and interesting deco- 
herence properties, which might provide a crucial advantage to using dot-confined hole spins for 
quantum information processing, as compared to electron spins. 
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I. INTRODUCTION 

The spin of a quantum-dot-confined electron is con- 
sidered a major candidate for the realization of solid- 
state-based quantum bits^ the basic building blocks for 
quantum information processing devices and, eventually, 
for a quantum computer £ One of the main obstacles to 
building these devices is the decay of spin coherence. The 
ultimate limit to the coherence time for electron spins in 
most quantum dots at low temperatures is set by the hy- 
perfine interaction with nuclei in the host material^i 4 *^ 
If no special effort is made to control this environment, 
the associated coherence times are quite short, typically 
on the order of nanoseconds 4^^ 

Very recently, several experiments have shown initial- 
ization and readout of single hole spins in self-assembled 
quantum dotsp ^ 10 : 11 and control over the number of 
holes in single gated quantunr dots,™ prerequisites for 
single-hole-spin dephasing-time measurements. Ensem- 
ble hole-spin dephasing times have recently been mea- 
sured in p-doped quantum wells. 13 Hole-spin coherence 
times in III-V semiconductor quantum dots are an- 
ticipated to be much longer than electron-spin coher- 
ence times due to a weak hyperfine coupling relative 
to conduction-band electrons! 14 i 15 i 16 i 17 i 18 i 19 i 20 ' 21 In the 
present work, we show that, in contrast, the coupling of 
a heavy hole (HH) to the nuclear spins in a quantum 
dot can be rather strong, potentially leading to coher- 
ence times that are comparable to those for electrons. 
However, in the quasi-two-dimensional (Q2D) limit, this 
interaction takes-on a simple Ising-like form, 



H = J2A h k sJ£, 



(1) 



where A\ is the coupling of the HH to the fc th nucleus, 



and If is the z- 



s z is the hole pseudospin-i operator, 
component of the k nuclear-spin operator 1^. The 
form of this effective Hamiltonian has profound conse- 
quences for the spin dynamics. Coherence times can be 



dramatically extended by preparing the slowly-varying 
nuclear field in a well-defined state ("narrowing" the 
field distribution) ^ i 22 ' 23 i 24 : 25 i 26 : 27 For an electron spin 
interacting with nuclei via the contact hyperfine inter- 
action, narrowing is effective only up to the time scale 
where slow internal nuclear-spin dynamics or transverse- 
coupling ("flip-flop") terms become relevant. Here we 
will show that heavy holes confined to two dinrensions 
have negligible flip-flops, potentially leading to signifi- 
cantly longer spin coherence times. The strong coupling 
of the HH to the nuclear spins is not due to confinement 
but is also present in bulk crystals, while the Ising-like 
interaction is a feature of Q2D systems. 

This paper is organized as follows: In Sec. [ill we write 
down the nuclear-spin interactions and derive an effective 
spin Hamiltonian for a quantum-dot-confined HH. In Sec. 
Mil we calculate the dynamics of the transverse HH spin 
for different external magnetic field directions. In Sec. 
IIVI we give estimates of the coupling strengths for the 
special case of a GaAs quantum dot. Conclusions and 
comparison to recent experiments can be found in Sec. 
IVl Technical details are deferred to Appendices [AUDI 



II. NUCLEAR-SPIN INTERACTIONS 



Hamiltonians 



For a relativistic electron in the electromagnetic field 
of a nucleus with non-zero spin at position Rfc, there are 
three terms that couple the electron spin and orbital an- 
gular momentum to the spin of the nucleus: the Fermi 
contact hyperfine interaction {h\ ), a dipole-dipole-like in- 
teraction (the anisotropic hyperfine interaction, h^), and 
the coupling of electron orbital angular momentum to the 
nuclear spin (/if). Setting H = 1, these interactions are 



described by the following Hamiltonians: 28 



Interactions in an atom 
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Here, 7s = 2/zb, 7^ = gj k [iN, Hb is the Bohr magneton, 
g Jfc is the nuclear g- factor of isotopic species jk, Hn is 
the nuclear magneton, r^ = r — R& is the electron-spin 
position operator relative to the nucleus, d ~ Z x 1.5 x 
10~ 15 m is a length of nuclear dimensions, Z is the charge 
of the nucleus, and rifc = r^/rfc. S and Lfc =rj.xp denote 
the spin and orbital angular-momentum operators of the 
electron, respectively. 

Nuclear-spin interactions are typically much weaker 
than the spin-orbit interaction. It is therefore appro- 
priate to form effective Hamiltonians with respect to a 
basis of eigenstates of the Coulomb and spin-orbit inter- 
actions. The 8x8 Kane Hamiltonian, which describes the 
band structure of a III-V semiconductor, provides such 
a basis * 29 ' 30 The Kane Hamiltonian is usually written in 
terms of conduction-band (CB) and valence-band (con- 
sisting of HH, light-hole (LH), and split-off sub-band) 
states. We derive an approximate basis of eigenstates in 
the HH sub-band by projecting the 8x8 Kane Hamilto- 
nian onto the two-dimensional HH subspace. 

To form effective Hamiltonians, we must approxi- 
mate the crystal-Hamiltonian eigenfunctions given by 
Bloch's theorem for a single band n: ^nkatr) = 

1 ^ikr„ 



/NX 



ikcr(r), where Na is the number of atomic sites 
in the crystal, and the Bloch amplitudes w„k<r (r) have the 
periodicity of the lattice. 

We will approximate the k = Bloch amplitudes 
UnOa-{ r ) within a primitive unit cell by a linear combi- 
nation of atomic orbitals (see Eq. (fj"2f . below). Near 
an atomic site, the CB Bloch amplitudes have approxi- 
mate s-symmetry (angular momentum 1 = 0), whereas 
the HH and LH Bloch amplitudes have approximate p- 
symmetry (I = 1). Adding spin, the ^-component of to- 
tal angular momentum of a HH is mj = ±3/2, whereas 
a LH has raj = ±1/2. In the Q2D limit, i.e., going 
from the bulk crystal to a quantum well (whose growth 
direction we take to be [001]), a splitting Alh develops 
between the HH and LH sub-bands at k = 0. We es- 
timate Alh — 100 meV for a quantum well of height 
a z ~ 5 nm in GaAs, much larger than the hyperfine cou- 
pling (see Appendix [A| . The splitting Alh is essential 
since it produces a well-defined two-level system in the 
HH sub-band, and we can restrict our considerations to 
the manifold of mj = ±3/2 states. 



Before addressing confinement in quantum dots, we 
illustrate that the interaction of an electron in a hydro- 
genic p orbital with the spin of the nucleus (which we 
choose to be at Rfc = 0) is generally non-zero. Moreover, 
when projected onto the manifold of mj = ±3/2 states, 
this interaction takes-on a simple Ising form. Although 
our final analysis will apply to any III-V semiconductor, 
we will take GaAs as a concrete example. 

The effective screened nuclear charges Z e ff "felt" by 
the valence electrons (in 4s and Ap orbitals) in Ga and As 
atoms have been calculated in Ref. [3l|. The 4s orbitals 
and Ap orbitals (with orbital angular momentum m^ = 
±1) are given in terms of hydrogenic eigenfunctions with 
the replacements Z — > Z c s by \&4oo(r) = Rio{r)Y§{Q, p) 
and ^4i±i(r) = Rn(r)Y 1 l (9, ip), respectively. Including 
spin and evaluating matrix elements of the Hamiltonians 
©-(SI with respect to hydrogenic 4s states leads to ef- 
fective spin Hamiltonians of the for m 28 ' 32 h\ s = A S S ■ Ifc 
and, due to the spherical symmetry of the wavefunction, 
/i| s = /i| s = 0. The same procedure with the Ap states 
leads to effective Hamiltonians h^ = (since p-states 
vanish at the origin) and h 2 p + h 3 p = A p s z I^, where 
s z = ±i corresponds to mj = ±3/2. A s and A p denote 
the coupling strengths of electrons in 4s and Ap orbitals, 
respectively. Evaluating all integrals exactly gives 



Ap = 1 f Z cS ( K ,Ap) 
A s 5 \Z cB (k,As) 



k = Ga, As. 



(5) 



Quite significantly, after inserting values of Z c ff from 
Ref. [3X| into Eq. {5]), we find that the ratio of cou- 
pling strengths is fairly large - on the order of 10%: 
A p /A s ~ 0.14 for Ga, and A p /A s ~ 0.11 for As. Since 
h\ and h\ do not contribute to the hyperfine interaction 
of an electron in an s orbital, research on hyperfine in- 
teraction for electrons in an s-type conduction band has 
focused on the contact term h\, neglecting the other in- 
teractions. The fact that h\, in contrast, gives no contri- 
bution for an electron in ap orbital has led to the claim 
that electrons in p orbitals (and holes) do not interact 
with nuclear spins. Eq. ([5]) shows that h 1 ^ and h\ can 
contribute significantly. Furthermore, while the interac- 
tion of a 4s-electron with the nuclear spin is of Heisen- 
berg type, within the manifold of mj = ±3/2 states, 
the interaction of a 4p-electron is Ising-like at leading or- 
der (virtual transitions via the mj = ±1/2 states may 
lead to non-Ising corrections). This result is a direct 
consequence of the Wigner-Eckart theorem and the fact 
that the Hamiltonians ©-(H]) can be written in the form 
hf = vf ■ Ifc , where v^ are vector operators in the electron 
(spin and orbital) Hilbert space. 



C. Interactions in a quantum dot 

We now return to the problem directly relevant to 
a HH in a two-dimensional quantum well. We con- 
sider an additional circular-symmetric parabolic confin- 
ing potential in the plane of the quantum well defin- 
ing a quantum dot. Neglecting hybridization with other 
bands, which we estimate to be typically on the or- 
der of 1% (see Appendix E}, the pseudospin states for 
a HH within the envelope-function approximation read 
|* ff ) = \4>\UYmo<j)\cr), where |</>;uhho<t) and \a) (a = ±) 
denote the orbital and spin states, respectively. The or- 
bital wavefunctions are given explicitly by (r|</>; uhhou) = 
Y / T?o0(r)wHHOcr(r). Here, vq is the volume occupied by a 
single atom (half the volume of a two-atom zinc-blende 
primitive unit cell), and (f>(r) = <fi z (z)4> p (p) is the enve- 
lope function. The radial ground-state envelope function 
is a Gaussian: 



w = -^i exp irw 



(6) 



where p = (x,y), p = \p\, I = l [1 + {B z /B f]- l l\ 
B z is the component of an externally applied magnetic 
field along the growth direction and B = & /ttIq where 
$0 = h/\e\ is a flux quantum. A typical dot Bohr radius 
of l = 30 nm gives B ~ 1.5 T. 

In a solid, the HH is delocalized over the lattice sites 
of the crystal. The nuclei do not interact solely with the 
fraction of the HH in the same primitive unit cell ('on- 
site' interaction), but also with density localized at more 
distant atomic sites (long-ranged interactions). We ne- 
glect the long-ranged interactions, which lead to correc- 
tions on the order of 1% relative to the on-site interaction 
(see Appendix [C| . If the envelope function varies slowly 
on the length scale of a primitive cell, we find (combining 
h\ and h\) A\ = A{ fc w o |0(R fc )| 2 , where 



A h - ijr ^sl Jk 



3 cos 2 fe + 1 

r 3 k (l + d/r k ) 



(7) 



p.c. 



Here, (• • • ) p . c . denotes the expectation value with respect 
to «HHo<r(r) over a primitive unit cell (a = ± give the 
same matrix elements), and 9 k is the polar angle of r^. 
The magnetic moment of a HH is inverted with respect 
to that for an electron. This results in a change of sign 
in Eqs. ©-(Bl and leads to the minus sign in Eq. ([7]) 
and of the values in columns (i) and (ii) of Table |TJ 
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FIG. 1: Dephasing of HH pseudospin states (solid lines in (a) 
and (c)). The decay is Gaussian for an out-of-plane magnetic 
field B z (a) (see Eq. (|10[) ) and given by a slow power law 
at long times (~ 1/Vi) for an in-plane magnetic field B x (c) 
(see Eq.|TT|). We have chosen |B| = 10 mT in (a) and (c). 
Magnetic-field dependences of the relevant coherence times 
are shown for a field that is out-of-plane (b) and in- plane (d). 
We have assumed g\\ = 0.04 (from Ref. |33|) and a zero-field 
lateral dot size lo = 30 nm and height a z = 5 nm, leading to 
TV = ttIqCIz/vo = 6.5 x 10 5 nuclei within the dot at B z = 0. 
We have taken vo = o|/8, where ah = 5.65 A is the GaAs 
lattice constant. The dashed lines in (a) and (c) show the 
dephasing of CB electron spin states in the high-field limit, 
|B| 3> (Je/g e ^B, where g e is the electron g-factor, and a e 
is obtained from Eq. <(9j by replacing A J h by the electron 
hyperfine coupling constants A\. 



between the mj = ±3/2 states via the LH sub-band 
are suppressed by ~ AJJ/Alh ^ 1- Hybridization with 
other bands can also lead to non-Ising corrections. For 
unstrained quantum dots, we find that these corrections 
are small: typically on the order of 1% of the values 
given in Table I (see also Appendix [AJ . Strain can lead 
to considerably stronger band mixing and, hence, to sig- 
nificantly larger non-Ising corrections to Eq. (jTJ). 



III. SPIN DECOHERENCE 

Now that the effective Hamiltonian (Eq. |T|)) has been 
established, we can analyze the dephasing of a HH pseu- 
dospin in the presence of a random nuclear environment. 
In an applied magnetic field, pseudospin dynamics of the 
HH are described by the Hamiltonian 



H = (b± 



(8) 



D. Non-Ising corrections 

There will be small corrections to the form of the ef- 
fective Hamiltonian given in Eq. UJ. Evaluating off- 
diagonal matrix elements of $5§ and Q with the ap- 
proximate Bloch amplitudes Q12[) yields non-Ising terms, 
whose associated coupling strengths Afc we find to be 
small: Aj^ < 0.06^. Higher-order virtual transitions 



where h z = J^k -^U-k ^ s ^ e nuclear field operator, b± = 
g±^BB z is the Zeeman splitting due to a magnetic field 
B z along the growth direction, and b» = gn/isBx is the 
Zeeman splitting due to an applied magnetic field B x in 
the plane of the quantum dot. g± and gi\ are the compo- 
nents of the HH g-tensor along the growth direction and 
in the plane of the quantum dot, respectively (we assume 
the in-plane g-tensor to be isotropic). 

If no special effort is made to control the nuclear field, 



the field value will be Gaussian-distributed in the limit 
of a large number of nuclear spins^ The variance for a 
random nuclear-spin distribution is (see Appendix [D]) 



(hi) 



1 

AN 



5>p'(p + i)(4) 2 



(9) 



where N = Trl 2 a z /vo is the number of nuclei within the 
quantum dot. The nuclear- field fluctuation a therefore 
inherits a magnetic- field dependence from I (see Eq. ([6])). 
A finite nuclear-field variance will result in a random 
distribution of precession frequencies experienced by the 
hole pseudospin, inducing pure dephasing (decay of the 
components of hole pseudospin transverse to B). 

First, we consider the case b\\ = 0. For a hole 
pseudospin initially oriented along the x direction, 
we find a Gaussian decay (see Fig. [1] (a)) of the 
transverse pseudospin in the rotating frame (s+) t = 
exp(-i6jj) ((s x ) t +i(s y ) t ): 



j 


Ai (neV) 


Ai (neV) 




(i) 


(ii) 


(iii) 


(iv) 


by Ga 
71 Ga 

75 As 


-7.1 
-9.0 

-8.2 


-13 

-17 
-12 


40 
51 
59 


74 
94 
89 


Aa = Yjj u i A a (a = e,h) 


-8.0 


-13 


52 


86 



TABLE I: Estimates of the coupling strengths for a HH 
(A 3 h ) and a CB electron (A J e ) for the three isotopes in GaAs. 
Columns (i) and (iii) show values obtained from a linear com- 
bination of hydrogenic eigenfunctions, using free-atom values 
of Z e g calculated in Ref . l3ll . Column (iv) gives the accepted 
values of Ai from Ref. 13a . Column (ii) shows the rescaled 
values from column (i) (see text). In the last row we give the 
average coupling constants weighted by the natural isotopic 
abundances: M39 Ga — 0.3, ^7i Ga = 0.2, za-s As = 0.5. 



halfway along the Ga-As bond by a linear combination 
of atomic orbitals, following Ref. I36l : 
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(10) 



This is the same Gaussian decay that occurs for 
electrons i^£ Here, since the magnetic field is taken to 
be out-of-plane, we must take account of the diamag- 
netic "squeezing" of the wavefunction. This squeezing 
affects the number TV of nuclear spins within the dot and 
hence, the finite-size fluctuation er. The coherence time 
T±(B Z ) = l/a{B z ) = rj_(0) [l + (B z /B ) 2 ]~ 1/4 then de- 
creases for large B z (see Fig. Q](b)). This undesirable ef- 
fect can be avoided for confined electron spins by generat- 
ing a large Zeeman splitting through an in-plane (rather 
than out-of-plane) magnetic field. This option may not 
be available for a HH where, typically, g\\ <C gj_. 

The situation changes drastically for an in-plane mag- 
netic field (b± =0). In this case, since the hyperfme 
fluctuations are purely transverse to the applied field di- 
rection, the decay is given by a slow power law at long 
times (see Fig. [1] (c)) and the relevant dephasing time 
increases as a function of the applied magnetic field (see 
Fig. [11(d)). In the limit 6|| ^S> a and for a HH pseudospin 
initially prepared along the i-direction, we find 



cos (but + \ arctan(t/r|| )) 
2[l + (-L)2U/4 



(11) 



The derivation of Eq. (fTTj) is directly analogous to that 
for the decay of driven Rabi oscillations in Ref. [34]. 



IV. ESTIMATES OF THE COUPLING 
STRENGTHS 

To estimate the size of A 3 h , we need an explicit ex- 
pression for the HH k = Bloch amplitudes. We ap- 
proximate UHHOo-(r) within a Wigner-Seitz cell centered 



UHHO<r(r) 



rGWS 



N a „<*„#£% (r + d/2 



+ v/l - ^K^(r - d/2] 



(12) 



Here, d = f(l, 1,1) is the Ga-As bond vector, a is 
the lattice constant, a v describes the relative electron 
sharing at the Ga and As sites in the HH sub-band, 
and N av is a normalization constant, chosen to enforce 
J ws d 3 r |uHHo<r(r)| 2 = 2, where the integration is per- 
formed over the Wigner-Seitz cell defined above. To 
simplify numerical integration, we replace the Wigner- 
Seitz cell by a sphere centered halfway along the Ga- 
As bond, with radius given by half the Ga-Ga nearest- 
neighbor distance. We find the electron sharing in the 
CB from the densities given in Ref. [35| to be a 2 ~ 1/2 
(see Appendix [BJ and assume the same (a 2 = 1/2) for 
Eq. (|12[) . Using Z e s for free atoms^i we evaluate Eq. 
([7]) with the ansatz, Eq. ([L2")) , by numerical integra- 
tion, giving the values shown in Table [U column (i). 
We check the validity of this procedure by writing the 
CB Bloch amplitudes as in Eq. (fT2)) . replacing the Ap- 
eigenf unctions by 4s-eigenfunctions. Evaluating the cou- 
pling constants for the CB (A J e ) from h\ gives the num- 
bers in column (iii). The accepted values of A 3 e from Ref. 
|35| are shown in column (iv) for comparison. Our method 
produces A\ to within a factor of two of the accepted val- 
ues. Our procedure, which relies on free-atom orbitals, 
most likely under- estimates the electron density near the 
atomic sites, which should be enhanced in a solid due 
to confinement. Assuming the relative change in density 
going from a free atom to a solid is the same for the CB 
and HH band, we rescale the results in column (i) by the 
ratio of the values in columns (iv) and (iii), giving column 
(ii). Due to the approximations involved, we expect the 
values in columns (i) and (ii) only to be valid to within 
a factor of two or three. 



V. CONCLUSIONS 

We have shown that the interaction of a quantum-dot- 
confined heavy hole with nuclear spins is stronger than 
previously anticipated. We have estimated the associated 
coupling strength to be on the order of 10 fieV in GaAs 
- only one order of magnitude less than the hypcrfinc 
coupling for electrons. However, the interaction turns 
out to be Ising-like which has profound consequences for 
hole-spin decoherence. Since no flip-flop terms occur in 
the effective Hamiltonian (Eq. ([T])), the main source of 
decoherence is given by the broad frequency distribution 
of the nuclear spins. Recent theoretical and experimen- 
tal studies have shown that state-narrowing techniques 
are capable of strongly suppressing this source of deco- 
herence, which makes the heavy hole an attractive spin- 
qubit candidate. 

Very recently, experimental results on hole-spin re- 
laxation in self-assembled quantum dots have been 
releasedi&ii Gerardot et al. report an extremely weak 
coupling of HH spin states, which is explained by our the- 
ory to be a direct consequence of the Ising-like nuclear- 
spin interaction (negligible flip-flop terms). Eble et al., 
in contrast, find very short hole-spin relaxation times on 
the order of 15 nanoseconds. This is due to the strong 
strain present in the particular dots used in this exper- 
iment, resulting in a considerable HH-LH mixing and a 
highly non-Ising interaction (large flip-flop terms). 
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APPENDIX A: HEAVY-HOLE STATES 

In this section we give details on our derivation of 
an approximate basis of heavy-hole (HH) eigenstates in 
a quantum dot. We will approximate the ground-state 
quantum-dot envelope function in the HH sub-band by 



(r) =<f> z (z)4> p (p), 
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where a z is the width of the confinement potential along 
the growth direction (for definition of the other symbols 
see Eq. (JBJ). We will then estimate the size of the split- 
ting Alh between the HH and the light-hole (LH) band 
and the degree of hybridization with the conduction band 
(CB), LH and split-off (SO) sub-bands. 



We start from the 8x8 Kane Hamiltonian given in 
Ref. l37| for bulk zinc-blende-type crystals, which is writ- 
ten in terms of the exact eigenstates (near k = 0) of an 
electron in the CB, HH, LH and SO band, usually de- 
noted by |l/2;±l/2) c , |3/2;±3/2)„, |3/2;±l/2)„, and 
|l/2;±l/2)„, respectively. We neglect terms that are 
more than two orders of magnitude smaller than the fun- 
damental band-gap energy E^ and perform the quasi- 
two-dimensional limit by assuming that a confinement 
potential has been applied along the growth direction. 
If the confinement potential is sufficiently strong (i.e., 
if the quantum well is sufficiently narrow), the energy- 
level spacing will be large and the electron will be in 
the ground state at low temperatures. Any operator act- 
ing on the ^-component of the electron envelope function 
can then be replaced by its expectation value with re- 
spect to the z-component of the ground-state envelope 
function. For the Kane Hamiltonian this means that we 
can replace powers of the z-component hk z of the crys- 
tal momentum Kk by expectation values. Assuming an 
infinite square-well potential of width a z confining the 
electron along the growth direction, the ground state is 
given by Eq. (|A3|) . Calculating the expectation value 
of k z and k z with respect to the ground state, we find 
(k z ) — and (k z ) = it 2 /a 2 . This allows us to write the 
Kane Hamiltonian in the following form: 
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where 
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Here, 



h 2 /2mo and toq is the free-electron mass, 



whereas m! is the effective mass of a CB electron. Fur- 



thermore, k± = k x ± ik y , 7' denote the Luttinger param- 
eters, P is the inter-band momentum, and Ago is the 
spin-orbit gap between the LH and the SO bands. Ex- 
perimental values for these parameters can be found in 
Table HI 

We assume a circular-symmetric parabolic confinement 
potential with frequency ujq in the xy-plane defining 
a quantum dot. Including a magnetic field along the 
growth direction, the ground state is approximately de- 
scribed by the Gaussian given in Eq. (|A2|) . The envelope 
function of the quantum dot is then the product of the 
in-plane and out-of-plane components, as given in Eq. 

In the quasi- two-dimensional limit, a gap 



h 2 i 2 



A LH = {B-C) = ^ fe + (k 2 y ) - 2(k 

mo V u 



(A5) 



develops between the HH and LH sub-bands, lifting the 
HH-LH degeneracy. Here, (• • • ) denotes the expectation 
value with respect to (|A1[) . The in- plane level spacing 
scales like ~ l/l 2 , where I is the dot Bohr radius. The 
in-plane level spacing is much smaller than the level spac- 
ing along the growth direction since, for typical dots, 
a 2 <C I 2 . Neglecting (k 2 ) and (k 2 ) compared to (k 2 ) in 

it 2 /a 2 for a square- well 



Eq. (|A5[) and inserting (k 2 ) 
potential, we estimate 

2n 2 i 2 h 2 



A 



LH 



a 2 m,Q 



100 meV 



(A6) 



for a z ~ 5nm, using 7 2 ~ 2.06 for GaAs (see Table |H| . 
The HH-LH splitting is thus much larger than the typi- 
cal energy scale associated with the hyperfine interaction 
(A e ~ 90peV for CB electrons in GaAs). 

To derive the approximate electron eigenfunctions in 
the HH sub-band of the quantum well, we start from the 
Kane Hamiltonian (|A4[) . We use quasi-degenerate per- 
turbation theory up to first order in l/£ (where £ stands 
for Eg, A L h, or A L h + A S o), taking H HH as the unper- 
turbed Hamiltonian. 37 This leads to a band-hybridized 
state of the form 



I* 



HH,hyb/ 



) = N<j } Xnl^na'jUnQa). 



(A7) 



Here, a = ±, (r|</w; u n0a ) = ^^(rju^W is the 
product of envelope function and k = Bloch amplitude 
in band n (CB, HH, LH or SO), the prefactors \° describe 
the degree of band hybridization, and J\f a enforces proper 
normalization. 

In first order quasi-degenerate perturbation theory, the 
hybridization with the CB and the LH and SO sub-bands 
is described by the interaction terms V\, V4, and V5 in 
Eq. (|A4p . respectively. We estimate the degree of hy- 
bridization by applying these operators to a two-spinor 
containing the in-plane ground-state envelope function 
(|A2[) of the HH sub-band. For the hybridization with 
the conduction band, we find (for B = 0) 



1 
E, 



■Vi 






A CB </>CBp 



\P) 
(P) 



(A8) 



P (eVA) 
E 9 (eV) 
Ago (eV) 



10.5* 
1.52* 
0.34* 



6.98 T 
2.06 f 
2.93 1 



TABLE II: Values of band parameters used in this section; 
*taken from Ref.—; f taken from Ref.—. 



where 

(t>CB P ±{p) = -=(V>io(p) ± iipoi(p))- (A9) 

Here, ipnm(p) = ipn(x)ip m (y) and tp n (x) is the n th 
harmonic-oscillator eigenstate. The envelope function 
of the admixed CB state is a superposition of excited 
harmonic-oscillator eigenfunctions. The prefactor 



A CB 



= ± 



P 



V2Egl 



(A10) 



determines the degree of sp-hybridization. Using values 
from Table[TT]and assuming a quantum dot with dot Bohr 
radius Iq ~ 30 nm (B = 0), we estimate A^ R ~ 10 -2 . 



Similarly, we estimate A LH ~ A so ~ 10 _ , assuming a 
dot height a z ~ 5 nm. The admixture of CB, LH and SO 
states to the HH state is thus on the order of 1% and has 
therefore been neglected in our considerations. 

We emphasize that sp-hybridization will lead to a cou- 
pling of the HH to the nuclear spins via the Fermi contact 
interaction @. Since the Fermi contact interaction is 
of Hcisenberg-type, sp-hybridization will directly lead to 
non-Ising corrections to the effective Hamiltonian given 
in Eq. (QJ. The size of these corrections is determined 
by the degree of sp-hybridization which is on the order 
of 1% (see above). 



APPENDIX B: ESTIMATE OF THE FERMI 
CONTACT INTERACTION 

In Eq. (fT2")h we have approximated the HH k = 
Bloch amplitudes within a Wigner-Seitz cell by a linear 
combination of atomic orbitals. Similarly, we approxi- 
mate the k = Bloch amplitude in the CB by 



wcBOo-(r) 



rGWS 



=N a la c ^ Q %(r + d/2) 



-V^#l$o(r-d/2) , (Bl) 



independent of a. Here, ^ 4 oo(r) = i?4o(r)y" °(6 l , ip), a c 
describes the relative electron sharing between the Ga 
and As atom in the Wigner-Seitz cell chosen to be cen- 
tered halfway along the Ga-As bond, and N ac normalizes 
the Bloch amplitude to two atoms in a primitive unit cell. 
The radial wavefunction depends implicitly on the effec- 
tive nuclear charges 2 , e ff(«, 4s), where k = Ga, As. 

We will estimate the relative electron sharing in the 
CB by calculating the electron densities at the sites of the 



nuclei from Eq. (IBID and comparing to accepted values 
taken from Ref. [35|. We will then estimate the Fermi con- 
tact interaction of a CB electron using free-atom effective 
nuclear charges taken from Ref. 31 (Z c ff(Ga, 4s) ~ 7.1, 
Z off (Ga,4p) ~ 6.2, Z off (As,4s) ~ 8.9, and Z off (As,4p) ~ 
7.4) and normalizing the Bloch amplitude over a Wigner- 
Seitz cell. 

We approximate the electron densities at the Ga and 
As sites within a primitive unit cell from Eq. (|B1[) : 



rfca = \u CBOa (-d/2)\ 2 ~ <,a»|*g£(0)| s 
d As = \u CBOa (+d/2)\ 2 ~Nl c (l 



o£)\*&> 



(0) 



(B2) 
(B3) 



We estimate the corrections to the right-hand sides to 
be on the order of 1% due to overlap terms. We take 
the ratio dGa/dAs and equate this with the ratio of the 
values from Ref. H d' Ga = 5.8 x 1CT 31 itT 3 and d' As = 
9.8 x 10~ 31 m~ 3 . This allows us to write a c as a function 
of the two effective nuclear charges: 



d' As f Z cS (Ga,4s) 

Z cff (As,4s) 



'Ga 



1-1/2 



(B4) 



Recalling that N ac normalizes the Bloch amplitude to 
two atoms over a Wigner-Seitz cell, we write 



N n 



\\ d3r 

- l Jws 


«c*S(r + d/2) 


9 1 


- \ 


/I - a^oo(r - d/2) 





1-1/2 



(B5) 



For all numerical integrations, we approximate the 
Wigner-Seitz cell by a sphere centered halfway along the 
Ga-As bond with radius equal to half the Ga-Ga nearest- 
neighbor distance. Inserting (|B4|) and (|B5|) into (|B2|) and 
(|B3|) . we solve the two coupled equations 



d Ga (Z off (Ga),Z off (As))- 
d As (Z cS (G&),Z eS (As)) 



d'c 

- d' t 



Ga 



As 



0. 
0, 



(B6) 

(B7) 



for the two effective nuclear charges. This yields 
Z c ff (Ga) ~ 9.8 and Z c s(As) ~ 11.0. Inserting these val- 
ues back into Eq. (|B4|) . we estimate the electron sharing 
within the primitive unit cell to be 



o; 



0.46. 



(B8) 



For comparison, inserting free-atom effective nuclear 
charges into Eq. (|B4[) yields a similar value: a' 2 ~ 0.54. 
Now we estimate the Fermi contact interaction of a 
CB electron starting from the free-atom effective nuclear 
charges Z e g(n, 4s) obtained from Ref. [3l|. We use a c ~ 
l/\/2 and normalize the k = Bloch amplitude to two 
atoms over a Wigner-Seitz cell, following Eq. (|B5|) . From 
the normalized Bloch amplitudes we estimate the Fermi 
contact hyperfine interaction by evaluating 



M 



2^0 



7s7jI m cbo<t(Rj)I' 



(B9) 



Here, R., = T-d/2 for Ga and As, respectively (j indexes 
the nuclear isotope). Evaluating for the isotopes in GaAs, 
this gives the values shown in column (iii) of Table Q] 

Replacing the Wigner-Seitz cell by a sphere with radius 
R s equal to half the Ga-Ga nearest-neighbor distance in 
our numerical integrations overestimates the expectation 
value of [r|(l + d/r fc )] _1 in Eq. 0. To estimate the 
error, we perform an integration over a sphere with radius 
R' s = (R s + i? max )/2, where R max denotes the radius of 
the smallest sphere that fully contains the Wigner-Seitz 
cell. From this, we estimate the relative error to be less 
than 30%. 



APPENDIX C: ESTIMATE OF THE 
LONG-RANGED INTERACTIONS 

In this section, we estimate the corrections to the 
HH coupling strength in Eq. J7|) due to long-ranged 
dipole-dipole interactions and long-ranged L • I interac- 
tions. To this end, we consider a single nucleus interact- 
ing with a HH that is delocalized over the lattice sites 
in the quantum dot. We start from the Hamiltonians 
given in Eqs. $3§ and ((3]). We define effective radii 
a e ff(«:,4p) = ao/Z c ff(K, Ap), where ao ~ 5.3 x 10~ n m 
is the Bohr radius. The effective radii define an approxi- 
mate length scale for the spread of the site-localized func- 
tions ^4 lo .(r) and are much smaller than the GaAs lattice 
constant a^ ~ 5.7 x 10 -10 m. The nucleus thus effectively 
'sees' sharp-peaked electron densities centered around 
the more distant lattice sites. We choose the nucleus 
to be at site R& and estimate the interaction with the 
electron density at more distant atomic sites by approx- 
imating the electron densities by 5- functions. Adding 
up contributions from h\ and /i|, we arrive at an effec- 
tive Hamiltonian describing the long-ranged interactions: 
Hfc = A^SzIfr, where Af r = J2i-i^k^h is tri e associated 
coupling strength and A™ — wo|^(R;)| 2 J^ 3 ^fc{^( r fc ~ 
~R-ki){h\ + h\)} describes the coupling of the electron den- 
sity at site R; to the nucleus at site R^ (R^ = Rj — R&)- 
In order to estimate the size of the long-ranged inter- 
actions relative to the on-site interactions, we take into 
account nearest-neighbor couplings for the long-ranged 
part, i.e., we replace J2i-i^k ~ * J2i=n n ■ ^ ne interac- 
tion with electron density located around more distant 
nuclei is suppressed by ~ ^/Rfi- Assuming that the 
quantum-dot envelope function varies slowly over the 
nearest- neighbor distance (</>(R;) — </>(Rfc) for I nearest 
neighbor of k) , we estimate the ratio of long-ranged and 
on-site interactions (Table I, column (i)) to be 



A h 



7 x 10" 3 , 



(CI) 



on the order of 1%, where Af r = ^4i r wo|^(RA:)| 2 . 

We remark that, in principle, the electron g-factor can 
deviate from the free-electron g-factor due to spin-orbit 
interaction. According to Ref. |39i . this renormalization 



is negligible for the on-site interaction, but could become 
relevant for the long-ranged interaction. However, for 
the estimate in Eq. (|Clj) . we have taken the free-electron 
g- factor. 



where 



/„ = Uo 2 5>(R fc )r 

k 



(D4) 



APPENDIX D: VARIANCE OF THE NUCLEAR 
FIELD 

Here we calculate the nuclear-field variance for a HH 
interacting with nuclei in a quantum dot. In particular, 
we evaluate 



(hi) 



(Dl) 



where (■ ■ • ) = TVj (p~j • • • ) indicates the expectation value 
with respect to the infinite-temperature thermal equilib- 
rium density matrix pi and we recall h z — ^ fe A\l^. For 
an uncorrelated and unpolarized nuclear state, we have 
VVw) = (4 2 > (#) = 0, k * k', which gives 



<r 2 =j2(A h k ) 2 mf 



(D2) 



Using ({ID 2 ) = P k (I jk +l)/3 for an infinite- 
temperature state, A\ = Ajf vo\4>(Rk)\ 2 , and assuming 
that the nuclear isotopic species with abundances Vj are 
distributed uniformly throughout the dot gives 



Assuming that the envelope function </>(r) varies slowly 
on the scale of the lattice, we replace the sum in Eq. (|D4|) 
by an integral: 

«o^|^(R fe )| 4 ->|d 3 r|0(r)| 4 . (D5) 



Inserting the envelope functions (|A2|) and (|A3|1 for a 
quantum dot with height a z and radius I and evaluating 
the integral in Eq. (|D5|) . we find 



P = 



3 1_ 

AN' 



(D6) 



Here, N is the number of nuclear spins within the quan- 
tum dot, given explicitly by 



N = 



irl 2 a z 



t'o 



(D7) 



a 2 = -I ^ Vj P(P + l)(Al) 2 , (D3) Inserting Eq. ||DBJ) into Eq. JD3]) directly gives Eq. ©. 
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